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Abstract 

We study the low-energy collective excitations of the neutron star inner crust, where a neutron 
superfluid coexists with a Coulomb lattice of nuclei. The dispersion relation of the modes is 
calculated systematically from a microscopic theory including neutron band structure effects. These 
effects are shown to lead to a strong mixing between the Bogoliubov- Anderson bosons of the neutron 
superfluid and the longitudinal crystal lattice phonons. In addition, the speed of the transverse 
shear mode is greatly reduced as a large fraction of superfluid neutrons are entrained by nuclei. 
Not only does the much smaller velocity of the transverse mode increase the specific heat of the 
inner crust, but it also decreases its electron thermal conductivity. These results may impact our 
interpretation of the thermal relaxation in accreting neutron stars. Due to strong mixing, the mean 
free path of the superfluid mode is found to be greatly reduced. Our results for the collective mode 
dispersion relations and their damping may also have implications for neutron star seismology. 
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I. INTRODUCTION 



The crust of a neutron star represents only about 10% of the star's radius and 1% of its 
mass but is expected to play a key role in various observed astrophysical phenomena such 
as pulsar glitches, quasi-periodic oscillations in Soft-Gamma Repeaters (SGR), and thermal 
relaxation in soft X-ray transients [1]. The outer crust is primarily composed of pressure 
ionized atoms arranged in a regular crystal lattice and embedded in a highly degenerate 
electron gas. With increasing density, electrons become relativistic and the rapid growth 
of their Fermi energy drives nuclei to become more and more neutron-rich due to electron 
captures 2|. Eventually, at a density ~ 4 x 10^^ g/cm^, some neutrons start to drip out 
of nuclei 3|, This defines the boundary between the outer and inner crust. "Dripped" 
neutrons in the inner crust are expected to become superfiuid below a critical temperature 
of the order of ~ 10^° K (see e.g. js]). Despite the absence of viscous drag, the neutron 
superfiuid can still be coupled to the crust due to non-dissipative entrainment effects arising 



a. 



As a matter of fact, it 



from elastic scattering of dripped neutrons by the crystal lattice 
has been shown that in some regions of the inner crust only a very small fraction of dripped 
neutrons participate in the superfiuid dynamics [8|, [Oj. Consequently, the vibrations of the 
crystal lattice are expected to be strongly coupled to the collective excitations of the neu- 



tron superfiuid 



10l-|l2j. Collective excitations are particularly important for understanding 



thermal and transport properties of accreting neutron stars with temperatures in the range 
T = 10^ - 10^ K (l3|. 

In this paper, we study low-energy collective modes with wavelengths large compared 
to the typical inter-nuclei distance. The existence of two longitudinal modes in the inner 
crust and the role of entrainment in determining the dispersion relations were first studied in 

n 

Ref. pj| using a hydrodynamic approach. In this pioneering study, long-range perturbations 
on the superfiuid flow induced by the lattice of nuclei were neglected. Here we show that they 
play a crucial role. The low-energy constants that depend on the microscopic properties of 
the inner crust are calculated in a consistent approach that properly incorporates the long- 
range correlations leading to entrainment effects, first discussed in Ref. [8|. We find that 
entrainment of superfiuid neutrons by crustal nuclei greatly reduces the velocity of the two 
transverse (shear) modes, and this in turn enhances their contribution to the low temperature 
specific heat of the inner crust. Entrainment effects also induce a strong mixing between 



2 



the longitudinal lattice phonons and the Bogoliubov- Anderson (BA) bosons |15l . Il6j of the 
neutron superfluid, splitting these modes into a high velocity global sound mode and a low 
velocity mode characterized by a relative motion between the neutron superfluid and the 
electron-ion plasma. These results should be also relevant for studies of global neutron star 
seismic modes with frequencies in the range 20 — 1000 Hz which could be excited during 
violent events such as giant flares in soft-gamma ray repeaters and binary neutron-star 
mergers. 

In the following section we deflne our notations and present order of magnitude esti- 
mates of the relevant length and momentum scales. In § IIIII we describe the low-energy, 
long-wavelength, collective excitation modes and their velocities. The § IIVI emploies the 



microscopic model of Ref. jo] to obtain quantitive values for these velocities. Damping of 
these modes is briefly considered in § |Vl In the § IVII we study how entrainment affects the 
inner crust speciflc heat, its thermal conductivity and its thermal relaxation time-scale. We 
flnally conclude in § IVIII 

II. BASIC NOTATIONS AND PHYSICAL SCALES 

In what follows, we will assume that the inner crust of a neutron star is a perfect crystal 
with a body-centered cubic lattice. Each crustal layer will be supposed to contain only one 
type of nuclide and will be characterized by Z the total average number of protons in the 
Wigner-Seitz (W-S) cell of the crystal lattice (a truncated octahedron), Aceii the total cell 
average number of nucleons, A the cell average number of nucleons bound inside nuclei and 
Ay, the cell average number of nucleons entrained by the solid crust. As shown in Ref. j^. 
Ay, is generally much larger than A and close to Acdi due to the Bragg scattering of unbound 
neutrons by the periodic potential of the crystal lattice, which manifests itself in neutron 
band structure effects. 

We will indicate hj n = Up+Un the total average baryon number density, which is the sum 
of the average proton density and average neutron density n„. The density of neutrons 
entrained by nuclei and the density of neutrons involved in superfluid motion will be noted 
as and respectively. As shown in Ref. sl, the density ra^ is generally much smaller 
than the density of "free" or "dripped" neutrons. Because of Galilean invariance, we 
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have 

nn = n^ + nl. (1) 
These densities are related to Z, A]^ and Acew by 

= n , (2) 

^cell 

and 



n. 



^cell 



(3) 



The ion number density ni is determined by 



"1 = 1—- (4) 

^cell 



As discussed in Refs. 
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12| . the mass density associated with lattice vibrations is given by 

pi = m {rip + n^) = A^mni, (5) 

where m is the nucleon mass (neglecting the small difference between neutron and proton 
masses), whereas the total mass density (neglecting the electron contribution) is 

p = mn = Aceii m rii . (6) 

The typical length-scale associated with the solid crust is the ion-sphere radius defined 



by 



A-nriiJ \ pi2 

where pi2 = p/(10^^g cm~^). The characteristic angular frequency and wave number of 
lattice vibrations is the ion angular plasma frequency 



OJr. 



'47r(Ze)2ni _ / 47re^nl 
Ay,m y {up + n^) m 



and the Debye wave number 

1/3 



g^ = (67rV)^/^^-^0.03f , ^'' ) im-\ (9) 
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respectively. The ion plasma temperature is defined by Tp = hup/k-Q (fee being the Boltz- 
mann constant). 

The ultra-relativistic electrons found in the inner crust of neutron stars with density 

n 

rig = Zrij, are almost uniformly distributed [17[ and are characterized by their Fermi wave 
number 

1/3 7 / 7 \ 1/3 / v/prn \ ^/^ 



I. 2 ,1/3 fzy ^7 fzy^^, fpi2z/5o_Y'\ 



(10) 



Small deviations of the electron distribution from uniformity are characterized by the elec- 
tron Thomas-Fermi screening wave number 



QTFe = \ UlTe^-— = \ kFe ^ 0.1 fcpe ~ — (H) 

" a/ie V TT ri \50/ 

where fig = fcpe is the electron chemical potential and a = /he ~ 1/137 the fine structure 
constant. 



III. LOW-ENERGY DYNAMICS OF THE NEUTRON-STAR INNER CRUST 



The equations governing the low-energy dynamics of a non-relativistic neutron superfluid 



.0 



12 



18| . The corresponding 



immersed in an elastic crust have been derived in Refs. 
normal modes of oscillation can be found by considering small perturbations of the densities 
and currents from their equilibrium values and solving the resulting linearized hydrodynamic 
equations. The first two of these equations arise from the conservation of neutron and proton 
numbers 

" + <V ■ + n\y ■ 5t;p = , (12) 



dt 

dSn 



^ + npV ■6vp = 0, (13) 



dt 

where 6nn and Sup are the perturbed neutron and proton densities respectively while 5Vn 
and 6vp are the perturbed neutron and proton velocities respectively. In the following we 
will consider oscillations characterized by wavevectors q <^ gxPe so that the crustal matter 
remains electrically neutral locally and Up = rif,. Treating the neutron-star crust as an 
isotropic solid and using i, j, k for coordinate space indices, the momentum conservation 



5 



can be expressed as 

™^n—^ + P^~^ + "-nVi^/in + LV iSfln - KV iUjj - 2SV j (^Uij - Sij^Ukk^ = , (14) 

where Uij is the strain tensor, (5/i„ the perturbed neutron chemical potential, S the shear 
modulus, K the bulk modulus of the electron-ion system 



2 I'dup ^ dfi. 



and the coefficient L given by 



rir. 



dup 



(16) 



takes into account the coupling of the neutron superfiuid to the strain field. The condition 
for neutron superfluidity is embedded in Josephson's equation 

^ + iv^,„ = 0. (17) 

ot m 

The normal modes have the form of plane waves that vary in space and time as exp[i(g • 
r — cut)], where q is the wave vector and u the angular frequency. In an isotropic medium, 
the normal modes may be separated into transverse and longitudinal ones. In the long 
wavelength limit g — t- 0, the normal modes all have a sound-like dispersion relation, with 
uj = vq, V being the mode speed. The speed of the two transverse lattice modes is given 
by HQ 

Vt = J-. (18) 

V Pi 

Due to interactions between neutron and proton densities and currents, the BA bosons of 
the neutron superfiuid with velocity are mixed with the longitudinal lattice phonons with 
velocity vi. Neglecting the coupling of the neutron superfiuid to the strain field, which 
is comparatively unimportant according to preliminary estimates, the resulting dispersion 



relation is given by H 



12| 



/ 2 2 2\f 2 2 2\ 2 2 2 

where the strength of the mixing is characterized by the parameter 



first introduced in Ref. 12|]. The velocity of tlie BA mode is 



m dun ' 

wfiereas tlie velocity of the longitudinal mode of the lattice is 



Pi 



(21) 



(22) 



In the neutron-star crust, the electron contribution to the bulk modulus dominates, and 
the ion contribution can be safely neglected (see e.g. Section 7.1 of 



approximately given by 
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QTFe 



Up 



rir, + m dUe 



ll). As a result, ve is 



(23) 



Solving Eq. (IT^ we find that the eigenmode velocities are given by 



V 



v± 



2^,2 



1 ± 



where 



2 

mix * 



(24) 



(25) 



The speed of the transverse lattice phonon in Eq. f fTSjl is unaffected by mixing and is 
approximately given by 20 1 



u„ ( Z 

vt ^ 0.4^ ^ 0.12 — 



1/3 



(26) 



Note that due to entrainment effects, the expressions (1?T]) . (125]) and for the veloc- 
ities of the BA bosons and lattice phonons differ from those obtained considering either a 
neutron superfluid alone or a pure solid crust, respectively. The self-consistent inclusion of 
entrainment is an important new element of this study. 

In the normal phase, any relative motion between the neutron liquid and the crust will 
be damped by collisions so that in the hydrodynamic regime ions, electrons and neutrons 
will be essentially co-moving. In this case, the Josephson's equation have to be replaced by 
the condition 6vn = Svp. As a result, only one longitudinal mode corresponding to ordinary 
hydrodynamic sound persists and its velocity is given by 



'K + AS/3 



(27) 



where K is the total bulk modulus of the crust. It is related to the bulk modulus K of the 
electron-ion system by 

2 ^l^n 



K = + 2n„L + < 



(28) 



Since S <^ K (see e.g. section 7 
as 



.lof fl 



), the sound velocity can be approximately written 



dp 



n m dn, 

The transverse mode velocity is given by 



Up Up dfXe _^ rin Un <9/i^ 



n m drir, 



(29) 



Vt 



(30) 



While the existence of two weakly damped longitudinal modes is unique to the superfluid 
phase, entrainment is fairly insensitive to superfluidity provided the pairing gap A -C /i„ |21| . 



which is the case in most of the inner crust 



22|-|24|. 



IV. MICROSCOPIC MODEL FOR THE INNER CRUST OF A NEUTRON STAR 

The evaluation of the velocities of the collective modes requires the knowledge of the 
susceptibilities defined by due/dfie and dun/dfin, and number densities Up and for each 
given baryon density n. Both Hp and dug/dfie are uniquely determined by the composition 
of the inner crust (i.e. the variation of the electron density with n) because electric charge 
neutrality requires Up = and since the electron gas is weakly interacting 

The composition of the inner crust, taken from Ref. jsl, was obtained from a non-relativistic 
Skyrme effective nuclear Hamiltonian solved using the fourth-order Extended Thomas- Fermi 
method with proton quantum shell effects added via the Strutinsky-Integral theorem. This 
ETFSI method is a high-speed approximation to the self-consistent Hartree-Fock equations. 
These calculations were carried out with the Skyrme force BSkl4 underlying the HFB-14 
atomic mass model [25!], which yields an excellent fit to essentially all experimental atomic 
mass data with a root mean square deviation of 0.73 MeV. In addition, this effective nucleon- 
nucleon interaction was constrained to reproduce the equation of state of neutron matter, as 
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TABLE I: Ground-state composition of the inner crust of a neutron star {Z, AceW, A as defined in 
Section |II]) , as obtained in Ref. ^ , for various baryon densities n /mass densities p. The effective 



number of bound nucleons was calculated including band structure effects in Ref. 



a. 



calculated by Friedman and Pandharipande 26[ using realistic two- and three-body forces. 



This equation of state is in good agreement with more recent ab intio calculations 27H29| 
at densities relevant to the neutron-star crusts. 

As discussed in detail in an accompanying paper jo], neutron band-structure calculations 
are needed to determine n'j^. Here, we note that the key ingredient is the single-particle (s.p.) 
dispersion relation Sak being the band index and k the Bloch wave vector) given by the 
solution of the Schrodinger equation with the periodic mean-field obtained self-consistently 
from the ETFSI method. The superfluid density was then found from the equation 



m 



247r3^2 



J2 [ Nkeak\dS^''\ 



(32) 



where dS^""^ is an infinitesimal area element of the piecewise Fermi surface associated with 
the a band. As described in Ref. js], in most regions of the inner crust only a small fraction 
of dripped neutrons contributes to the superfluid density due to Bragg scattering so that 
<^ or equivalently ~ Aceii- Note that unbound (bound) neutrons with density ra^ 
(respectively — n^) are characterized by s.p. energies Eak lying above (respectively below) 
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TABLE II: Properties of collective modes in the inner crust of a neutron star. The velocities (u^, 
vi^ Cs, Vt, V-, and f+) and the mixing parameter gmix ai^e defined in Section HIH c is the speed of 
light. Values in parenthesis are obtained by neglecting entrainment. The ratios of the mean free 
path of the longitudinal modes to that of the unmixed longitudinal lattice phonon are shown in 
the last two columns. 

the largest value of the periodic mean-field potential. Results are summarized in Table [H 

The neutron chemical potential is determined by the neutron band-structure from the 
equation 



deD{e) 



where D{e) is the density of neutron s.p. states defined by 



(33) 



(34) 



where the fc-space integration is taken over the first Brillouin zone. Differentiating Eq. fl33|) 
with respect to Hn thus yields the neutron susceptibility 

drir 



, dD(e) 
de- 



(35) 



Because nuclei in the inner crust are neutron-saturated, the neutron susceptibility is essen- 
tially independent of the neutron bound states except possibly in a small region close to 
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neutron drip. For the reasons explained in Ref. 30|, the density D{e) of neutron unbound 
states in a given region of the inner crust is well approximated by the density of s.p. states in 
uniform neutron matter for the corresponding density of dripped neutrons. Using these 
approximations, the velocity of the BA mode in the inner crust can be expressed as 




(36) 



where f ^ is the velocity of the BA mode in pure neutron matter at the density associated 
with the crustal layer under consideration. This latter velocity is given by 31 1 



3 



(37) 



where vp is the Fermi velocity in pure neutron matter at the density while Fq and 
Fi are the corresponding dimensionless Landau parameters, whose expressions for Skyrme 



interactions can be found in Ref. 



32| . We have evaluated w^, using the same Skyrme effective 



interaction BSkl4 as that used to determine the equilibrium composition of the crust. 
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FIG. 1: (Color online) Speeds (in units of c) of the longitudinal (left panel) and transverse (right 
panel) collective excitations in the inner crust of a neutron star. Dotted curves show results with 
neither mixing nor entrainment, dashed curves include effects due to entrainment only and solid 
curves include in addition the effects due to mixing. 

The speeds of the collective modes in the inner crust of a neutron star are shown in Fig. [H 
and listed in the Table [Tll Entrainment modifies the spectrum, v,/,, ve and Vt are all signifi- 
cantly reduced (compare dotted and dashed curves), and mixing leads to a strong splitting 
between the longitudinal eigenmodes (note the difference between speeds of the lowest and 
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highest eigenmodes). With increasing density, a strong suppression of the plasma frequency 
due to entrainment leads to rapid decrease in the velocity of transverse and longitudinal 
lattice phonon modes. Mixing between longitudinal modes leads to high velocity eigenmode 
with velocity f+ and low velocity mode with velocity v^. The v_ mode is predominantly 
the superfluid phonon (BA) mode near neutron drip and transforms to a mode with a large 
lattice component at the crust-core boundary. The mode with velocity v+ is a pure lattice 
mode at neutron drip and transforms to being a mode which is predominantly a superfluid 
mode at the crust-core interface. 



V. DISSIPATION 



Lattice phonons couple strongly to electrons and easily excite electron-hole pairs in the 



dense electron gas. This Landau damping of lattice phonons has been studied in Ref. [33 1 
and an approximate result of the lattice phonon mean free path was obtained. The mean 
free of a thermal phonon that contributes to thermal conductivity was found to be 

Gvr 1 FiTjT) ^^^^^ ( AQ\^'' ( F{T,/T) 

where 



Aiph = ^ 72.5 - ' r,eii , (38) 

Ze^ V Aph-e \Z J \V Aph_e 



Aph-e = In (-^ - (l- and 7 = qD/kpe , (40) 



.7/ 2 V 4 

and V is average velocity of the lattice phonon. Note that for simplicity we have neglected 
corrections due to the Debye- Waller factor and the nuclear form factor to the Coulomb 
logarithm Aph_c. Such corrections tend to increase the mean-free path and Eq. f l38|) must 
therefore be viewed lower limit. 

Our interest here is to investigate the mean free path of the superfluid phonon mode 



in the inner crust. In Ref. 34] it was shown that phonon-phonon and phonon- impurity 
scattering were negligible compared to the dissipation that arose due to mixing with the 
lattice phonon. The superfluid phonon mean free path, without the inclusion of entrainment 
effects, was found to be much larger than that of the lattice phonons because mixing due to 
the density interaction was weak. In the following, we include effects due to entrainment, 
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which is now known to be the dominant contribution to the mixing parameter gmix, and 
we show that the mean free path of the superfluid mode is greatly reduced due to strong 
mixing. Incorporating this into the dispersion relation in Eq. (IT^ we obtain 

{u' - vl q'){uj' - 2ir, u - v} q') = gl^ , (41) 

where = Vi/Xg and Xe is the mean free path of the lattice phonon in the limit of weak 
damping (F^ ^ viq). In general, Xi ^ Aiph as the latter is an average mean free path 
more closely related to the mean free path of the transverse thermal phonon. Nonetheless 
it provides an order of magnitude estimate. 

Mode mixing induces an indirect coupling between the superfluid BA bosons and elec- 
trons. Because the longitudinal modes contain an admixture of superfluid and lattice 
phonons, the damping of lattice vibrations due to electron-hole excitations naturally leads 
to a finite damping of both modes. In a small region of the crust in the vicinity of the 
neutron-drip transition where (/^j^, "^^l — v^^ the modes are not strongly mixed, and using 
the fact that vi ^ f^, we can obtain from Eq. fHTj) the analytic relation 

between the mean free paths of the superfluid and lattice modes. In other regions mixing is 
strong and damping associated with each eigenmodes is found by solving Eq. f HT]) . Although 
an analytic solution exists, it is cumbersome to write down explicitly. We present numerical 
values for the ratio of the mean free paths A+/A^ and X^jXi where X± are the mean free 
paths of the eigenmodes in the last two columns of Tabled It is meaningful to calculate these 
ratios without specifying Xi because it is independent of in the weak damping limit. From 
the table we see that the mean free path of the mode with a large superfluid component is 
large near neutron drip, but decreases rapidly due to mixing when vn ~ v^. In the bulk of the 
inner crust both modes have comparable mean free paths and this behavior is qualitatively 
different from that observed in Ref. 3J] where entrainment was neglected and mixing was 



found to be weak except in a narrow region close to resonance. 



VI. IMPLICATIONS 



X-ray observations of accreting neutron stars in low-mass X-ray binaries have recently 
proved to be very useful for probing neutron-star interiors. The accretion of matter onto 
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the surface of the neutron star triggers thermonuclear fusion reactions. Under certain cir- 
cumstances, these reactions can become explosive, giving rise to X-ray bursts and super- 
bursts 35j . The ignition conditions of these thermonuclear flashes depend sensitively on the 



thermal properties of the crust. Valuable information on neutron star crusts can also be ob- 
tained from the thermal X-ray emission in quiescence following a long outburst of accretion 
during which the crust has been driven out of its thermal equilibrium with the core 36 |. 
The thermal relaxation between the accreting and quiescent stages has been monitored for 
the four qu asi-p ersistent soft X-ray transients KS 1731—260 37(|, MXB 1659—29 



J1701-462 
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XTE 



and EXO 0748-676 



40| . Numerical simulations of these phenomena have 



shown that the cooling is very sensitive to the properties of the neutron-star crust [36 
In particular, the thermal relaxation time of the crust is approximately given by 
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r ~ {ARY 1 



Rc^ J 



Cv 

K 



(43) 



where AR is the crust thickness, R is the radius, and M is gravitational mass of the neutron 
star, while Cy and k, are the average heat capacity and thermal conductivity in the density 
range between ~ O.lficc and n^c where Ucc = 0.08 fm~^ is crust-core transition density. 
The thermal relaxation of hot newly-born neutron stars could also shed light on the crust 
properties. However, such very young neutron stars have not been observed yet, being 
presumably obscured by their expanding supernova envelope. 

The inner crust heat capacity is the sum of contributions from the quasiparticle excita- 
tions of the electron gas and neutron liquid, and from the collective excitations described 
above. In what follows we describe these contributions to the volumetric crustal heat capac- 
ity. At densities above ~ 10^ g cm~^, electrons can be treated as an ideal relativistic Fermi 
gas so that their heat capacity is simply given by (/eb^ ^ /Ue) 

1 f^l 



CI 



rknT. 



(44) 



3 {hcY 

The heat capacity of non-superfiuid degenerate neutrons (for k-^T <^ /i„) is similarly given 

by 

1 



(45) 



where D{fin) is well approximated by the^ density of states in uniform neutron matter at 
the density n^, i.e. -D(/i„) ~ rnkpnl {t^^Y [30,]. This neutron contribution is enormous and 
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FIG. 2: (Color online) Left panel: Debye temperatures 0/) of the longitudinal collective excitations. 
Dotted curve show values with neither mixing nor entrainment, dashed curve include effects du to 
entrainment only and solid curves include in addition the effects due to mixing. Right panel: Debye 
temperature 0£i of the transverse collective modes and ion plasma temperature Tp = {fi/k-Q)^p- 
In both panels, the light blue band delimits the range of T below which nuclei crystalize (using 
Vc = 180 to 220). 

will always dominate in the layers where neutrons are normal. Once superfluidity sets in, 
however, Cy is strongly suppresssed and becomes negligible when the temperature is much 
lower than the critical temperature T" 23|, ■ Given the density dependence of the neutron 
■So gap there are „,rly two regies, just above the neutron^drip pomt arrd possibly i,r the 
deepest part of the crust, where Cy is relevant (see, e.g., [13]). 

The heat capacity associated with a collective excitation having a dispersion relation of 
the form u = vq is given by 



Srir. 



dx- 



X e 



(46) 



[e-^ - 1)^ 

with xn = Qd/T, Qd = (^/^b)q'd'^ being the Debye temperature of the collective mode. 
At low temperatures, T ^ 6d such that ^ 1, one has the standard Debye result 

3 A 4 / 1 rri \ 3 



"YColl 



271^ fknT 



An' 



rii 



(47) 



15 \ hv J 5 \hqDV J 

while at high temperatures, T ^ Qo when ^ 1, one obtains the classical result C™" = 
rii. At low enough temperatures (T ^ 6d and T ^ T") the heat capacity of the crust is. 
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FIG. 3: (Color online) Heat capacity of electrons (e) and transverse lattice phonons (tph) and 
longitudinal excitations (— and +) in the inner crust of neutron stars at T = 10'' K, with mixing 
and entrainment effects (solid lines) and without (dotted lines). In the absence of mixing, the 
longitudinal modes are the Bogoliubov-Anderson superfluid phonons (sph) and the longitudinal 
lattice phonons (Iph). For comparison, is also shown the normal neutron contribution (n), but it is 
strongly suppressed by superfluidity except in the shallowest and densest parts of the inner crust 
where the neutron ^Sq pairing gap becomes vanishingly small. 

hence, approximately given by 



We plot in Fig. |2] the Debye temperatures of the four collective modes and the various 
contributions to Cy are displayed in Figs. I3]l4l andS] for three typical temperatures of as- 
trophysical interest. Fig. [2] shows that these three temperatures are well below the Debye 
temperature of the f _|_ mixed mode while at T = 10^ K the f _ and the two degenerate trans- 
verse modes are in the classical regime but at T = 10^ and 10'' K all modes are well into 
the Debye regime. This suggests that entrainment and mixing will not affect the thermal 
relaxation of newly-born isolated neutron stars but could be important for accreting neutron 
stars. 
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FIG. 5: (Color online) Same as Fig. [3] for T = 10^ K.. 



Figs. [3l m anciS] show that the heat capacity of non-superfluid neutrons would largely 
dominate over all collective modes, but becomes insignificant once neutron superfiuidity 
sets in, i.e., in most of the inner crust. Overall, the transverse lattice mode contribution Cy 
to the heat capacity, dominates at T = 10^ K and 10^ K, while electrons dominate at 10'' K 
due to the linear temperature dependence of Cy compared to the dependence for Cy. Cy 
in not affected by entrainment at T = 10^ K, since the transverse modes are in the classical 
regime, but at T = 10^ K and 10^ K it is increased by almost one order of magnitude in most 
of the crust. Notice that at K without entrainment Cy would be comparable to Cy while 
it clearly dominates once entrainment is take into account. Moreover, the heat capacity 
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of the longitudinal mode is increased by several orders of magnitude by entrainment, and 
mixing. In particular, the contribution of the lowest mixed mode becomes even comparable 
with Cy at high temperatures. 

Because entrainment modifies the spectrum of collective excitations, it also impacts the 
heat transport in the crust. The thermal conductivity is generally governed by electrons. 
Changes of phonon velocities alter the electron-phonon process hence also the electron ther- 
mal conductivity. The conductivity is mainly limited by the Umklapp process in which an 
electron simultaneously Bragg scatters off the lattice, and emits a transverse phonon {45, 46 |. 
Since the scattering rate scales as where v is the phonon velocity, and vi ^ Vt, processes 
involving longitudinal phonons are typically negligible. This observation also permits us to 
reliably estimate the changes in the electron mean free path due to entrainment. First, we 
note that the electron-phonon scattering rate depends on the electron Fermi momentum fcpe, 
the ion plasma frequency Up and Vt (see [isl of a discussion). Since Vt oc uip/qD, it follows 
that effects due to entrainment on the scattering are entirely incorporated through its ef- 
fects on Up. It therefore suffices to employ existing fitting formula developed in earlier work, 
but with a suitably reduced value of Up due to entrainment. In Fig. [6] we plot the electron 
thermal conductivity with and without entrainment effects included. As anticipated, the 
conductivity decreases with entrainment simply reflecting the fact that it is easier to excite 
lower velocity transverse phonon modes. 

Although electrons dominate heat conduction under normal conditions, phonons can 
contribute either at high temperature when C y^^ > Cy or when large magnetic fields suppress 



electron conduction transverse to the field 



33 



3^. In t 



le inner crust^the lattice and 



superfiuid phonons contributions were estimated in Ref. [33| and Ref. [3J], respectively. 
From kinetic theory and in the case where phonon conduction is diffusive (rather than 
convective) the thermal conductivity is given by 

3 

where C™'' is the heat capacity, v is the velocity and A is the mean free path of each collective 
mode. Entrainment alters the thermal conductivity through these three factors. The larger 
specific heat associated with lower velocity transverse modes implies that their contribution 
to the heat conduction is proportionately enhanced. In addition, since oc the smaller 
Vt acts to further increase the conductivity, and the combined effect is to increase the earlier 
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FIG. 6: Electron thermal conductivity in the crust with (solid) and without (dashed) entrainment 
effects included, at five densities n = 0.0003, 0.001, 0.01, 0.02, and 0.08 fm~'^ as labeled on the 
curves. The minimum of /tg occurs at T ~ O.lTp, marked by a dot on the corresponding curve: 
Kg oc in the quantum regime, T <^ Tp, while it only weakly increases with T in the classical 
regime at higher temperatures. 



estimate of Ref. j33| by the factor {Ak/Ay/'^. 

The effects on the superfluid phonon contribution is more complex because mixing is 
strong throughout the inner crust except in the vicinity of neutron drip. It is only mean- 
ingful to discuss heat diffusion due to eigenmodes, and in general there are two competing 
effects due entrainment. At first entrainment lowers the velocity of the mode with a larger 
superfluid component and increases its heat capacity, but with increasing density this in- 
crease is overcome by strong mixing which dramatically reduces the mean free path. Since 
A+ and A_ are of the same magnitude as A^, and because t>+ ^ Vt and f _ ^ Vt their con- 
tribution to heat transport is typically negligible. This new result implies that superfluid 
modes may play a smaller role in heat transport in magnetars than anticipated in Ref. 34| . 
and it is likely that the enhanced heat conduction due to the transverse mode will dominate 
in much of the inner crust. 



VII. CONCLUSIONS 

A large fraction of dripped neutrons in the inner crust of a neutron star are entrained by 
nuclei and move with them, owing to elastic (Bragg) scattering of neutrons by the crystal 
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FIG. 7: Ratio of the heat capacity Cy to the electron thermal conductivity Kg at five temperatures, 
as labeled. Cy includes the ion as well as the electron contributions, the neutron part is neglected. 
Continuous lines show values when entrainment, through its modification of Tp, is taken into 
account, while in the values for dotted lines it is neglected. 

lattice . This non-dissipative entrainment induces a strong coupling between the su- 
perfiuid and lattice dynamics, and is shown to impact the spectrum of low-energy collective 
excitations of the inner crust. Superfiuid and longitudinal lattice phonons are found to be 
very strongly mixed, and the speed of transverse lattice modes is greatly reduced thus lead- 
ing to a significant enhancement of the crustal specific heat at temperatures above ~ 10^ K. 
This, combined with entrainment induced reduction in the electron mean free path, entails 
an increase of the heat diffusion time in the crust, especially for temperatures in the range 
10^ — 10*^ K encountered in quasi-persistent soft X-ray transients. This warrants the need to 
take into account entrainment effects in the interpretation of the observed thermal relaxation 
in these accreting neutron stars. 

Shear modes in neutron-star crusts with velocity in the range Vt — 10~^ — 10~^ c have 
been proposed to play a role in the interpretation of quasi-periodic oscillations (QPOs) 
observed in giant flares from soft-gamma ray repeaters 47|]. The fundamental frequency of 
the global shear mode is given by VLq ~ Vtjl'nR^ R being the neutron-star radius, and Vt an 
appropriate average of the shear velocity in the inner crust, where the mode energy mainly 
resides 48|. Since entrainment lowers Vt by a factor of about 2 — 3 in most of the inner 



crust, our results suggest that Qq is too small to account for the observed QPO frequencies 
in the giant flares 49|]. It is also likely that the existence of the low- velocity longitudinal 
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eigenmode in the coupled superfluid-solid inner crust may be relevant to interpret global 
oscillation modes. 

However, there are several issues that deserve further attention before one can draw 
quantitative conclusions from our study. The possible presence of nuclear "pastas" in the 
deep regions of the inner crust [Jj], which has been neglected here, would reduce the effects 

6| and change the temperature dependence of the specific heat at 



of Bragg scatterin 
low temperatures 



50[ | due to the low dimensionality of these configurations. Besides, the 



composition and the properties of neutron-star crusts may differ from those of cold- catalyzed 
matter that we have considered in this work. We anticipate that quantum and thermal 
fluctuations of nuclei about their equilibrium positions, crystal defects, impurities and more 
generally any source of disorder would presumably reduce the number of entrained neutrons. 
Quantitative estimates of all these effects is beyond the scope of this work. 
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